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1  Summary 

A  mathematical  equation  has  been  developed  to  calculate  the  open-circuit  voltage  for  a  direct 
methanol  fuel  cell  in  the  absence  of  electrode  poisoning  and  methanol  crossover  by  the 
application  of  thermodynamic  fundamentals  of  the  phase  and  electrochemical  reaction  equilibra. 
The  equation  accounts  for  the  effect  of  nonideal  behavior  of  fluid  phases  in  the  anode  and 
cathode  compartments  of  a  fuel  cell  on  the  open-circuit  voltage.  The  developed  equation  also 
accounts  for  the  effect  of  different  total  pressures  in  the  anode  and  cathode  plenums  on  the  cell 
voltage.  For  the  ideal  behavior  of  the  fluid  phases  in  the  anode  and  cathode  compartments  the 
equation  normally  reduces  to  the  Nemst  equation.  The  developed  equation  is  capable  of 
predicting  the  open-circuit  voltage  at  any  desired  fuel  cell  temperature. 

2  Introduction 

Modeling  activity  in  conjunction  with  the  experimental  work  is  deemed  to  economically 
aid  in  the  development  of  an  optimal  design  of  a  direct  methanol  fuel  cell  (DMFC)  for  aerospace 
and  ground  transportation  vehicle  applications.  To  this  end,  research/development  activity  by 
the  principal  investigator  (P.I.),  Sarwan  S.  Sandhu  is  given  below. 

1.  The  principal  investigator  (P.I.)  developed  a  theoretical  scheme  to  predict  the 
reversible  cell  voltage  in  the  absence  of  electrode  poisoning  and  methanol  crossover  through  a 
polymer  electrolyte  membrane  of  a  DMFC,  for  example,  Nation  "  perfluorosulfonic  acid  polymer 
and  polybenzoxazole  or  polybenzimidazole  based  electrolyte  membranes.  The  scheme  based  on 
the  application  of  thermodynamic  fundamentals  of  the  phase  and  electrochemical  reaction 
equilibria  predicts  the  reversible  voltage  at  any  fuel  cell  temperature.  The  effect  of  nonideal 
behavior  of  fluid  phases  in  the  anodic  and  cathodic  fluid  compartments  of  a  fuel  cell  on  the 
reversible  voltage  is  accounted  for.  The  simulation  is  capable  of  predicting  the  reversible  cell 
voltage  as  function  of  reactant  feed  composition  and  different  total  pressures  in  the  anode  and 
cathode  plenums.  The  developed  reversible  voltage  equation  reduces  to  the  Nernst  type  equation 
for  the  ideal  fluid  phase  behavior  of  the  fuel  and  oxidant  feeds.  The  predicted  reversible  cell 
voltage  can  be  used  as  an  ideal  standard  to  which  the  actual  voltage  of  a  direct  methanol  fuel  cell 
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operating  at  a  given  set  of  reactant  feed  composition,  temperature,  and  pressure  conditions  can  be 
compared  to  evaluate  the  fuel  cell  performance. 

A  senior  year  undergraduate  student,  Mr.  R.  Owen  Crowther,  typed  the  entire  theoretical 
development  by  the  P.I.  for  the  reversible  cell  voltage  prediction.  The  P.I.  guided  Owen  to 
develop  a  computer  code  to  simulate  the  entire  set  of  mathematical  equations  to  generate  the 
reversible  cell  voltage  data  as  function  of  temperature,  fuel  feed  composition,  and  anode  and 
cathode  side  total  pressures. 

The  actual  theoretical  development,  the  computer  code,  and  the  code  generated  numerical 
data  in  the  form  of  tables  and  plots  are  presented  in  Report  #1  entitled:  “Fuel  Cell  Project:  Direct 
Methanol  Fuel  Cell:  Theoretical  Formulation  of  Reversible,  Open-Circuit  Voltage  for  a  Direct 
Methanol  Fuel  Cell”  by  Sarwan  S.  Sandhu,  Department  of  Chemical  &  Materials  Engineering, 
The  University  of  Dayton,  Dayton,  Ohio  45469-0246. 

2.  The  P.I.  developed  the  first-attempt  mathematical  equations  to  predict  the  ionic 
current  density  and  the  species  (H^CFfiOH,  FLO)  molar  fluxes  for  their  transport  through  a  solid 
polymer  electrolyte  or  polymer-ceramic  material  composite  membrane  of  a  DMFC  at  the  steady 
state,  isothennal  conditions  by  the  application  of  fundamentals  of  transport  phenomena.  The 
developed  transport  equations  are  presented  in  three  sets.  In  the  transport  set  I,  the  equations  are 
based  on  the  extension  of  moderately  dilute  solution  theory.  These  equations  involve  the  Fick’s 
law  type  mass  diffusivities  that  are  known  to  be  strongly  dependent  on  the  species 
concentrations.  The  transport  equation  sets  II  and  III  are  valid  for  very  strong  solutions.  The 
transport  equation  set  II  is  based  on  the  application  of  the  generalized  Maxwell-Stefan  equations. 
The  transport  equation  set  III  is  based  on  the  application  of  the  Lars  Onsager’s  irreversible 
thennodynamic  approach  to  transport  processes.  The  developed  equations  account  for  the  effect 
of  the  voltage,  pressure  and  the  species  concentration  gradients  on  their  transport  fluxes.  These 
equations,  describing  species  transport  fluxes  through  a  solid  polymer  electrolyte  membrane  are 
to  be  coupled  with  the  equations,  yet  to  be  developed,  describing  mass  transfer  through  the 
electrode  porous  backing  layers  and  species  mass  transport  with  electrochemical/chemical 
kinetics  in  the  porous  electrode  layers.  The  coupled  equations,  describing  the  various 
phenomena  occurring  in  a  DMFC  in  operation,  would  be  of  immense  significance  in  the 
evaluation  of  performance  and  design/development  of  a  DMFC.  The  usefulness  of  the 
developed  transport  equations  may  also  be  seen  in  that  they  can  be  employed  to  experimentally 
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determine  the  transport  properties  of  a  membrane  such  as  permeability  of  the  fluid  mixture  and 
effective  mass  diffusivities  of  the  species  by  designing  appropriate  fuel  cell  experiments. 

The  developed  transport  equations  by  the  P.I.  were  typed  by  Mr.  Sarath  Krishnan,  a 
graduate  student  in  the  chemical  engineering  program  at  the  University  of  Dayton.  Very 
recently,  Mr.  R.  Owen  Crowther  (who  is  now  a  graduate  student  in  our  chemical  engineering 
program)  has  been  guided  by  the  P.I.  to  simulate  the  first  set  of  transport  equations  in  the  form  of 
a  computer  code  to  generate  numerical  data  on  current  density  and  species  (H+,  CH3OH,  H20) 
molar  fluxes  through  a  polymer  electrolyte  membrane  of  a  DMFC.  He  is  expected  to  start 
developing  the  computer  code  in  the  near  future. 

The  actual  theoretical  development  of  the  transport  equations  is  presented  in  Report  #2 
entitled:  “Fuel  Cell  Project:  Direct  Methanol  Fuel  Cell:  Theoretical  Formulation  of  Transport 
Fluxes  of  species  (H+,  CH3OH,  H20)  through  a  Solid  Polymer  Electrolyte  Membrane  (PEM) 
of  a  Direct  Methanol  Fuel  Cell,”  by  Sarwan  S.  Sandhu,  Department  of  Chemical  &  Materials 
Engineering,  The  University  of  Dayton,  Dayton,  Ohio  45469-0246. 

3.  The  P.I.  developed  a  set  of  first-attempt  equations  to  predict  the  reversible  as  well  as 
irreversible  or  actual  power  production  from  a  continuously  fed  DMFC  operating  at  steady-state 
conditions.  This  development  is  based  on  the  fundamentals  of  classical  thermodynamics  and  the 
species  and  overall  material  balances.  The  scheme  requires  the  information  on  the  reactant  feed 
stream  composition,  temperature,  and  pressure  conditions  at  the  inlets  of  a  DMFC,  fractional 
conversion  of  methanol  via  electrochemical  process  for  the  production  of  electric  power,  its 
chemical  oxidation  at  the  catalytic  surface  of  the  cathode,  and  the  amounts  of  methanol,  water, 
and  hydrogen  ions  transported  through  the  solid  polymer  electrolyte  membrane  from  the  anode 
side  to  the  cathode  side  of  a  DMFC.  The  development  accounts  for  the  nonideal  behavior  of  the 
fluid  phases. 

This  development  is  not  in  the  typed  form  at  the  time  of  reporting.  Mr.  Sarath  Krishnan, 
a  graduate  student  in  the  chemical  engineering  program  at  the  University  of  Dayton  has  gladly 
accepted  to  understand  the  entire  scheme  of  equations  and  to  present  them  in  the  fonn  of  a  typed 
report.  Also,  he  would  develop  a  computer  code  under  my  direction  to  generate  numerical  data 
on  the  reversible  and  actual  electric  power  production  as  function  of  input  parameter  values. 

4.  The  P.I.  studied  and  carried  out  partial  analysis  of  the  open  literature  on  the 
research/development  of  direct  methanol  fuel  cells  depending  on  availability  of  time  for  this 
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equally  important  activity.  Based  on  the  information  acquired  from  the  literature  on  the  direct 
methanol  fuel  cells  and  his  knowledge  of  the  basic  and  engineering  sciences  relevant  to  the 
research/development  of  a  DMFC,  the  P.I.  developed  an  outline  of  a  comprehensive  research/ 
development  program  to  develop  an  efficient  DFMC  energy  conversion  system.  The  outline  has 
been  entitled,  “Direct  Methanol  Fuel  Cell  (DMFC)  Research/Development  -  Theoretical  and 
Experimental  Integrated  Approach”.  A  copy  of  this  outline  was  provided  to  Dr.  J.P.  Fellner 
of  the  Power  Division  of  the  WPAFB,  Ohio. 

5.  The  effort  by  Sarwan  S.  Sandhu  (P.I.),  his  students  and  Dr.  J.  P.  Fellner  on  the 
research/development  activity  has  resulted  in  the  following  presentations  and  proceedings 
publications. 

(a)  “Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling:  Open-Circuit  Voltage 
Equation”  in  the  Proceedings  of  the  8th  International  Conference  on  Electrical  and  Electronic 
Products,  Vol.  34  (2002),  pp.  233-253,  Product  Safety  Corporation,  Sissonville,  West  Virginia; 
The  Conference  held  at  the  Greenbrier,  White  Sulphur  Springs,  West  Virginia;  January  14-16, 
2002. 

(b)  “Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling:  Reversible  Open  Circuit 
Voltage  Equation  and  Species  Flux  Equations”,  presented  at  the  8th  International  Symposium  on 
Polymer  Electrolytes”  held  at  the  Eldorado  Hotel,  Santa  Fe,  New  Mexico,  May  19-24,  2002. 

(c)  “Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling,”  presented  at  the  2002 
AIChE  North  Central  Regional  Student  Conference.  “A  Climate  Change,”  Chem-E-Car 
Competition,  University  of  Michigan,  February  7-9,  2002. 

3  Methods,  Assumptions,  and  Procedures 

Standard  chemical  engineering  thermodynamic  principles  were  used  to  develop  the  reversible 
open-circuit  voltage  for  a  direct  methanol  fuel  cell.  Open-circuit  voltage  equations  based  on  gas 
and  solution  thennodynamics  for  ideal  and  non-ideal  cases  were  derived  and  their  effect  upon 
open-circuit  voltage  was  determined. 
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4  Results  and  Discussion 


The  results  and  discussion  from  Report  #1  entitled  “Theoretical  Formulation  of  Reversible  Open- 
Circuit  Voltage  for  a  Direct  Methanol  Fuel  Cell”  are  given  below. 

4.1  Development  of  Performance/Design  Equations  for  a  Direct 
Methanol  Fuel  Cell 

4.1.1  Equilibrium  Open-Circuit  Voltage  Full  Equation 


Figure  1.-  Schematic  Representation  of  a  Direct  Methanol  Fuel  Cell  (not  to  scale). 

Legend 

x  =  distance  measured  from  the  anode  toward  the  cathode 
AP  =  Anode  plenum 
CP  =  Cathode  plenum 

A(P)=  A  porous,  platinum-catalyzed  electrode  (e.g.  Pt  or  Pt/Ru  alloy),  anode  labeled 
as  p  phase  (solid),  laminated  to  the  membrane  (M) 
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C(P')=  A  porous,  platinum-catalyzed  electrode  (e.g.  Pt  or  Pt/Ru  alloy),  cathode 
labeled  as  [V  phase  (solid),  laminated  to  the  membrane  (M) 
a,a'=  Platinum  (solid)  lead  wires 

M=  Membrane  phase  functioning  as  electrolyte  medium  and  separator  between  the 
electrodes.  This  membrane  phase  may  be  consisted  of  a  polymer  electrolyte  such 
as  Dupont  Nafion  perfluorosulfonic  acid  polymer  supported  by  porous 
polybenzoxazole  (PBO).  Or,  it  may  be  composed  of  a  ceramic  electrolyte 
in  composite  with  a  polymer  such  as  PVDF  (polyvinylidene  fluoride) 
to  enhance  the  ionic  conduction  by  providing  additional  transport  pathways 
between  the  ceramic  electrolyte  and  the  polymer  material  boundaries. 


Here  fuel  is  methanol  (CH3OH)  present  as  mixture  with  water  (liquid  phase)  at  the 


concentration  level  Canon. ap 


g  -  mold  If )l  I 
liter  of  solution 


or  rncmoH,  ap 


g  -  mold  If )I  l 


kg  of  solvent  H 20 

in  the  anode  plenum.  It  is  assumed  that  the  methanol  concentration  in  the  porous  anode  is 
uniform.  Also,  oxidant  species  concentration  and  product,  water,  species  concentration  are 
assumed  uniform  in  the  porous,  cathode  electrode.  Furthermore,  it  is  assumed  that  the 
isothermal  (i.e.  constant  temperature  T)  conditions  prevail  in  the  system  with  the  total  pressures 
in  the  anode  and  cathode  plenums  at  P)  and  Pc‘,  respectively. 

The  cell  reactions  are: 

At  anode: 

(1) 

phase) 


CH3OH  (1)  +  H20  (1)  <->  C02  (g)  +  6H  (1)  +6e 

(in  the  liquid 


phase  at  x=xA) 


At  cathode: 


3/2  02  (g)  +  6H+(1)  +6e  <->  3H20  (1) 


(in  the  liquid  (in  p '  phase) 
phase  at  x=xc) 


The  overall  cell  reaction: 


(2) 
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(3) 


CHj,OH  (1)  +  3/2  02(g)  <->  C02(g)  +  2H20  (1) 

I  I 

(in  the  anode)  (in  the  cathode)  (in  the  anode)  (in  the  cathode) 

The  approach  given  in  the  Electrochemical  Systems  [  J.S.  Newman;  Prentice  Hall  (1991 
edition)]  is  employed  here.  The  phase  and  electrochemical  reaction  equilibria  are  used.  Say,  the 
anode  and  cathode  electric  potentials  are  Oa  and  Oa  ,  respectively. 

The  difference  between  the  electron  electrochemical  potentials  at  the  cathode  and  anode 

is 

juae'  - juae_  =  (juf™™1’0'  +  Ze  FOa')-(^ceh_emical’a  +  ze_FOa)  (4) 

Because  jUcJ,em,cal’a'  =  ^m,Cai,a  por  qlc  identical  platinum  lead  wires;  therefore,  Eq.  (4)  becomes 
Me-  -  Me-  =  Ze-F^a  ~  ® “ )  (because  ze  =  - 1 )  (5) 

=  -F(0“'-0“)  (5a) 

=  -FU 

Therefore,  FU  =  -  <Da )  =  Me-  -  Me-  (6) 

where  U=  open-circuit  (or  reversible)  cell  potential 

F=  Faraday’s  constant  =96  487  Coulomb  per  equivalent=  magnitude  of  the  charge 
on  lg-mole  of  electrons  or  protons  (H+). 

Mae_  =  Me  (phase  equilibrium  with  respect  to  electron  as 

an  electrochemical  species)  (7) 


At  anode,  the  condition  is  of  the  reaction  equilibrium 

A  A  A  s'  A  s  p 

MCH30H  +  MhIO  =  Mcoi  ®Mh+  +  vMe- 
Therefore, 


m:  = 


f  A  A  _  A  A 

Mchioh  +  Mh20  Mco2 


A 

■Mm 


From  Eqs.  (7)  and  (8a),  the  following  expression  is  obtained. 


(8) 


(8a) 
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(8b) 


where  are  lbc  chemical  potentials  of  CH3OH,  H2O,  CO2  in  their  respective 

phases  at  the  anode;  and  /j^+  is  the  electrochemical  potential  of  H+  in  the  liquid  phase  in  the 
porous  electrode. 

Now,  jLi'f  =  fif  (phase  equilibrium)  (9) 


At  the  cathode,  the  condition  of  the  reaction  equilibrium  is 


3 

A02  +  ^Mh+  +  ~  3/JH20 


c  3  c  c 
3 Mh2o  ~  A02  ~~  ^/uH+ 

- V - 


(10) 


Therefore, 


B'  / 1  C  1  c  \  C 

Fe-  =  (t  A//20  _  T  A02)  ~  + 

2  4  | 

(11) 

Using  this  result  in  Eq.  (9),  the  following  equation  is  obtained. 

a'  _  /  1  C  1  C  \  C 

/4~  —  V^FhIO  ~^Fo2>  Fh+ 

(12) 

Substitute  for  //"  and  //“  from  Eqs.  (8b)  and  (12),  respectively,  into  Eq  (6)  to  obtain  Eq.  (13): 

Vi  1  A 


FU  =  F[$>a'  -®a)-- 


(a  a  a  \ 

McHWH  +  Mmo  —  Fc02 


J 


~Fh+ 


1  ,c  -I  ,c 

2  Fh  20  4  Fo 2 


-Mh+  (13) 
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FU  =  F(oa  -  <£")  = 


(a  A  i 

Mchwh  Mco2  ,  ,.c 

6  6  Fm 


f  ..A 


+ 


.C  \ 


Mu  20  Mh  20 


+  (/£-<)  04) 


J 


( MACH20H  .  1  ..c  ^ 

6  4^°2 


'+'  +(Mh+  ~  Mh+)  (15a) 


McQ2  ,  Mh2Q(1)  Mh2QQ)  ,  (,.C 

H*  '~7/4 


FU  =  F(oa'  -  r)  =  - 
v  ;  6 

F[/  =  f(o“'  -<£“)=- 


A  3^  C 

Mchioh  +  ~ZMo2 

2  y 


^  (/^C02  Mh20(1 )  +  3 Mh20(1 )  )+  ( 


y  3  c  a  ,  , 

McH30H(l)  +  ~ZMo2  —  \Mc02  +  \TMh20(1)  ~  Mh2o 
V  2  7 


A  ^  'b"c  )!+kc.-^.; 

(16) 


(assuming  methanol  and  water  in  the  liquid  phases  at  the  electrodes) 

[  chemical  potential  of  a  species  in  a  gas  phase  mixture  at  T,  P]  is  expressed  as: 


jui{  )  =  //;)  (@  T,  P°  =1  bar-standard  state  pressure)  +  RT  In 


\pC  J 


(17) 


where  fugacity  of  species  i  in  the  gas  mixture  at  (T,Pl)  with  composition  y,’s  (mole  fractions)  is 
given  by 

f,=biPt  (Ha) 

where  (j)i  =  (fugacity  coefficient  of  the  component  i  in  the  gas  mixture  to  account  for 

the  nonideal  behavior  of  the  species  i  in  the  gas  mixture  if  the  gas  pressure  is  greater  than 
10  bar). 


Substituting  the  above  result  into  Eq.  (17)  leads  to: 


Mi(g)  =  MaS)  (@  T,  P°[standard  state  pressure=lbar]  +  RT ln| 


h- 

TlSl  yyO 

V  r  J 


(17b) 


Here,  juf  (@  T,  P°[=lbar])  =  chemical  potential  of  pure  species  i  at  the  temperature  T 

and  pressure  P°=lbar  where  it  follows  ideal  gas  behavior; 
ij)t  =  <f):  (T,  Pl,  y,)=  fugacity  coefficient  of  component  i  in  the  gas  phase  mixture  to 
account  for  the  nonideal  behavior; 
v,  =  mole  fraction  of  component  i  in  the  gas  mixture. 
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jUi(l)  [chemical  potential  of  a  charge  neutral  species  in  a  liquid  phase  mixture  (T,P,Xj  s)], 
following  [J.S.  Newman,  Electrochemical  Systems  (1990  edition],  is  given  by 


Mm  =  hrinl,,,,]-  [»rlnX,J+  RT  In*,,, 


=  RT\nA°(l)+RT\n 


f ;  '\ 
\Am  ) 


-  jui(l)  \of  purei  at  T,P'  \+  RT\nar  i  (18) 

where  ar  i  =  relative  activity  of  species  i  in  the  liquid  -  phase  mixture; 

Xi(l)  =  absolute  activity  of  the  component  species  i  in  the  liquid-phase  mixture; 

A°U)  =  absolute  activity  of  the  pure  species  i  in  the  liquid  phase  at  the  liquid  phase 
mixture  temperature  and  pressure. 

According  to  Smith  et  al.  [Introduction  to  Chemical  Engineering  Thermodynamics,  McGraw- 
Hil,  Inc.,  (1996  edition)]. 

pm  =  pf  (@T,P°  =  1  bar)  +  RT\na,  (19) 

where  activity  of  species  i  in  the  solution  is  given  by 

'vt(pt  ~ P° ) 


a,  =  YiX , 


V  fi  J 


=  (ylxl)exp 


RT 


(19a) 


where  yt  activity  coefficient  of  species  i  in  the  liquid  phase  mixture;  it  is  a 


function  of  temperature,  pressure,  and  composition  (xjs), 

f  =  fugacity  of  pure  i  in  the  liquid  phase  at  the  mixture  temperature  and  pressure 
conditions, 

f°  =  fugacity  of  pure  i  in  the  liquid  phase  at  the  mixture  temperature  but  pressure; 
P°(=  1  bar), 

Vi  =  molar  volume  of  pure  i  in  the  liquid  phase  at  temperature  T. 


Combining  Eqs.  (19)  and  (19a)  leads  to: 

Man  =  k /)  i@T’  P°  =  1  bar)  +  V/P'  -  P°)\+  RTHytxt) 

(19b) 

=  p\ °m\ofpure  i  at  T,P\  +  RT  ln(;/ixi ) 

(19c) 

On  comparison  of  Eqs.  (18)  and  (19c),  one  finds 
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ar.  =  yixi  (19d) 

Equation  (19b)  is  used  for  methanol  and  water  in  the  liquid  phase  mixture  to  obtain: 


Mch3oh(i )  —  Mch3oh(i)  (@  —  lbar) 

+  VcH,OH(n(P' -  P°) 


(20a) 


+  RT  \n(yl 


CH,OH(l)  ACH3OH(l)) 


Mh,o(1)  =  Mh,o(1)(@T>P°  =lbar) 


+  vHiom{P,-P°) 

+  RT\n(yHiOU) 

xH2o(iy 


(20b) 


where  VCH  OH(l) ,  VH  G(/)  are  the  molar  volumes  of  pure  CH3OH  and  H20  in  the  liquid  phase, 
respectively,  at  temperature  T. 

By  the  use  of  Eq.  (17b),  the  chemical  potentials  of  oxygen  and  carbon  dioxide  are  given  as: 


Ma,(g)  =  Mo,<g)  (@R  p°  =  1  bar)  +  RTln 


0o2yo2  po 


Mcoyg)  =  Pcoxg)  (@T,P°  =  lbar)  +  RTln 


J 

5 

r  CO,  y  C02  uo 
V  r  J 


(20c) 

(20d) 


Therefore,  |/^c//3o//(o  +  ^  P°2  j  —  Pen 3o//(/)(@T,T>  —  lbar)  +  VcH  }OH(l)  (P>A-  P°) 


+  RTln(;/( 


A  1 
CH  3OH(l)ACH3OH(l)> 


+  -<  (@T,P°  =lbar)  +  RTln| 


j'.cA 


iC  c 
Vo  Jo. 


I uco2  +  (3 


C  _  A  j!  _  A 

Mh2oo)  Mh20(1)J\  -  Pco 


)}=//Aft  (@T,  P°  =  1  bar)+  RTln 


2  ^  u2  no 

V  r  J 


f  p',A  \ 


lA  A 

rcojco 


2*  7”)  O 


(20e) 


+  3//  *  (@  T,P°=  1  bar)  +  3  VH,(/)  (Pl’c  -  P°  )  +  3RTln(y  "  ^  ) 


c 


H20{1) ' 


h2o(/) 


H20(/)  H  O(|)  ■ 


-Vh2O(0(rM  -P°)-RTln(rA0(/)x^rtm) 


H20(/)'  H2OQ) ' 
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=  Mco2  +2<0(/)  +  VH20(/)[(3P'-c -P'a)-2V°] 


+  RT  In 


1A  A 

rco ,  y  co 


fLU;'  po 

~A  A 

7  H20(l)XH20(l) 


Substituting  the  results  from  Eqs.  (20e)  and  (20f)  into  Eq.  (16)  to  obtain: 
FU  =  f(ok  -Ok) 


Mch3oh{1) (@T)  +  - // *  (@ T )  +  VCHjOH(!)  (P'’A-P°) 


+  f?rln<!  (/CvOH<l)XCH,OH0))(.<^O2yO2  ) 


j)  -Aco, 


-f?r  In 


pt,A 

2A  ,A  r 

Vco.yco,  po 

A  A 

Yh20(1)XH20(1) 


(20g) 


Fu  =  F  r  -a>° 
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f-1 

V  6y 


A  _3_  A  _  A  _  9  A 

McHjOhu)  +  ~^"o2  Mco2  ZMh20(1) 


{''c,o„m(p'-A  -p°)-y«,om^p'-c  -p,‘  -2P°)} 


+  RT  In 


Ych 


3OH  (l)XCH  3OH  (/) 


>  t,c 


,  C  C 

ro,  y  o 


2  y  u2  T)  ° 


I  A  A 

r  CO  2  y  CO 


f  p  t  ,A  \ 


v  P°  j 


Y  h2o{1)x  h2o(1) 


,  ,  c  ,4  \ 

+  l/v-  -/v; 


FU  =  F  ®K  -® 


(<£x  -®x) 


(20h) 


v6y 


{  A  _3  A  _  //A  _9  /A  I 

Mch3oh(I)  +^Po2  L1co2  zMh2o(1)( 

+  L«,r-f)-^,pi>'c-pu  -zf”)} 


+  i?r  In 


(fen 


A  A  A 

0//(/)/  H2O(/)XCH3O//(/)XH20(/) , 


3 

^  ,  P'-C  ^2 

2c  c  ■* 

'02  y  02  J-yQ 

v _ _  y 

A 

iA  A  1 

<Pco2yCo2  p0 


(21) 


N0W’  <+„ 


^electrochemical  potential  of  H+  in  the  polymer  elcctroytc  ' 
at  the  cathode  -  electrolyte  interface 
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(22a) 

condition  of  the  phase  equilibrium  with  regard  to  H 1  ions; 
note  that  H+  ion  is  the  mobile  ion  in  the  polymer 
electrolyte 

^electrochemical  potential  of  H+  in  the  polymer  electro lyte^ 
v  at  the  anode  -  electrolyte  interface  y 

(22b) 

Insert  for  pcH+(l)  and  ju^+(l)  from  Eqs.  (22a)  and  (22b)  into  Eq.  (21): 

FU  =  f(<Dx'  -Ok) 


A  A-E-I 

Mir  =  /V 


v6y 


A  _i__a_a_9a 

I  McH,OH(l)  +  ~2^02  Rc02  ZRh20(1) 


+  I  v, 


(PuA  -po^_vH  (3Pr’c  -Pl’A  -2P°) 


t,C  r)t,A 


CH3OH(l ) 


h2o(1) 


+  RT  In 


(/ CH3OH(l )Y  H20U)XCH30H(1)XH20(1)  f  <t>o2yco, 


u  2  no 


pt,A  \ 

1a  a  ■* 
rco,yco2  po 
V  r  J 


ps„-r-,p-\ 


(23) 


Following  the  approach  given  by  J.S.  Newman  (Electrochemical  Systems,  p.  83,  Prentice  Hall 
(1991)),  the  electric  potential  (also  called  quasi-electrostatic  potential)  in  the  solution  of  the 
membrane  phase  is  defined  as: 

juH+  =RT\ncH+  +zh+FO  (24a) 

Then, 
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PC^EI  -vt*-1 


RT In  +z  +jf(«Dc  £  7 -O'4  £  0 

nA-E-l  H*  \  / 


=  RT\w  P/f,  +  z  +F( ®c  E  1  -  ®A  E  ') 

„  A-E-I  H  V  / 


=  ^ln  +  V 


Q>AEI) 


(24b) 


where  mH+  =  molality  of  H  in  the  solution  with  water  as  solvent, 
c,,+  =  molar  concentration  of  H4  in  the  solution, 

P soiv  =  solvent  (water)  density, 

V  =  +1’ 

O  =  electric  potential. 

From  Eq.  (24b),  the  expression  for  (/z^' '  -  )  is  inserted  into  Eq.  (23)  to  obtain 

FU  =  Fjr'-o*) 


{  A  .  A  A  r\  A  I 

McH3OH(l)  +  ~  ^c°i  ~ZrlH20(l)  I 


!  +{rc„,oWt(p'A -n-VmoVP1* -p‘  A  -2F”)) 


tt}  ~t  I  q  ^  ^ 

+  K1  1m  y CH30H^Y h2o(1)xch3oh(1)x  h2o(1) 


/C  c  ± 

0o2yO2 


jiA  A  * _ 

'02yo2  po 


mc +E  1 

+  RT^^hrr  +ZH 

J 


+  zh+f(®c~e-1  -Qae/) 
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4.1. 1.1  Simplification  of  the  Full  Equation 

Equation  (25)  shows  that  the  cell  potential  is  dependent,  for  the  selected  standard  states  of  the 
species  (CH3OH,  O2,  CO2,  H20)  at  the  cell  temperature  T,  on 

•  total  pressures  at  the  electrodes,  the  standard  state  pressure  (P°); 

•  pure  molar  volumes  of  CH3OH(l),  FFO(l)  in  the  liquid  phase; 

•  mole  fractions  of  CH3OH,  H2O  in  the  liquid  phase  mixture  at  the  anode  and  the  mole 
fractions  of  02(g),  ycQi ,  and  C02  (g),  y  Aco  ,  at  the  cathode  and  anode,  respectively; 

•  activity  coefficients  of  CFEOFIQ,  Ych.ohu)  ,  ar|d  H20(1),  Yh,ou)  at  the  anode; 

•  fugacity  coefficients  of  02(g),(f/o  ,  and  C02(g),(J)A0  ,  in  the  gas  phase  at  the  cathode 
and  anode,  respectively; 

•  the  hydrogen  ion  concentrations,  mc~R~\  mA~E~' ,  in  the  membrane  at  the  cathode  and 

anode  interfaces  and  on  the  electric  potential  difference  (o  (  /_/  —  O  1  7/)  i  n  the 

electrolyte  membrane  between  its  interfaces  with  cathode  and  anode. 

Simplifications  of  Eq.  (25): 

(1)  If  the  potential  variation  across  the  membrane  is  ignored,  especially  for  the  case  of  no 
poison  reaction  such  as  the  parasitic  oxidation  of  methanol  at  the  cathode,  i.e. 

(<f>c  E  1  -  cp'4~£~y)  =  0  in  Eq.  (25),  the  open-circuit  potential  should  then  be  given 


by 


fu  =  f(«dcc-occ) 


6 


|  A  _l_A_A_9A 

Mco, 


+ 


{vc„,oh^(PLA  -P°)~  P,£  -  P'  ‘  -  2 P°)\ 


+  RTln 


(fcH: 


OH(l)Y H20(1)XCH,0H(1)XH20(1)  . 


-|3/2  M 


1C  c 

to  Jo. 


2  '  w2  O0 


■ >t,A 


lA  A 
TCO-y  y  CO 


2  ^  L  U2  T)0 


(  C-E-I  \ 

m 


+  RT  In 


m 


H+ _ 

A-E-I 


\  H+  J 


(26) 
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Comment:  Determination  of  the  cell  potential  using  this  equation  requires  information 

C-E-7 

m 

on  — in  addition  to  other  parameters.  These  can  be  detennined  if  one  has  the 
m 

H+ 

thennodynamic  information  of  the  partition  or  distribution  coefficient  defined  as 

M 

Kh+  = - ^ -  (26a) 

m/r,  (in  aqueous  solution) 

(obtained  from  the  phase  equilibrium  data) 


In  case  one  cannot  evaluate  the  last  tenn  in  Eq.  (26),  one  may  assume  the  last 
tenn  in  this  equation  to  be  negligible.  This  is  equivalent  to  the  assumption  that  Vm  +  in  the 
membrane  =  0  under  the  cell  open-circuit  conditions.  This  equation  is  then  reduced  to: 


FU  =  F  ®x  -® 


(®x  -®K) 


{  A  3  a  _  A  _0  .A  I 

McHiOmi)  +  Mco2  z"H2oo)f 


]_ 

6 


+ 


-P°)-VH,m(3P'c-P'-A-2P°)} 


+  RTln 


f  r 


[yL 


A  A  A 

OH(l)/  H20(l)XCH,0H(\)XH20U) , 


t,C 


2c  c  P 
Kyo,  po 


3/2  V 


"it,  A 


I A  A 

rco.y CO 


2  r>o 


(27) 


•  The  activity  and  fugacity  coefficients  of  the  species  can  be  detennined  using 
the  methods  given  in  the  standard  thermodynamic  books: 

Smith,  et  ah,  Introduction  to  Chemical  Engineering  Thermodynamics,  Chapters 


10-14;  The  McGraw-Hill  Companies,  Inc.,  5th  Edition  (1996); 

Tester,  J.W.,  and  M.  Modell,  Thermodynamics  and  Its  Applications,  Chapters  9,  1 1,  12, 
15;  Prentice  Hall;  1997  edition. 

J.S.  Newman,  Electrochemical  Systems,  Chapter  4  (for  ionic  activity  coefficients), 
Prentice  Hall  (1991  edition). 

•  Standard  state  chemical  potential  data  at  a  reference  temperature  and  methods  to 
obtain  them  at  any  other  desired  cell  temperature  are  also  provided  in  the  standard 
thennodynamic  books.  Also,  the  other  property  values  such  as  V;  [pure 
component  molar  volume]  can  be  obtained  from  the  references  given  above. 


For  the  case  of  dilute  solution  of  CH30H(1)  in  H20(l) 


r  a  , mo1 

e-§A  CCH3OH(l)  <0-1 


V 


and 


J 


total  anode  and  cathode  pressures  being  <  10  bars,  one  may  assume  that  the  liquid 
and  gas  phases  to  follow  ideal  solution  and  ideal  gas  behavior,  respectively.  With  this 
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assumption,  rw3omn  =  Yh2ou)  =  1;  K  =  kj2  =  1  equation  (27)  becomes 


FU  = 


|  A  3  a  _  A  _  9  ,A 

\Mch3ohu)  +  2  Ao2(g)  >“co2(g)  zMh2o(1) 


+  {vc„j„m,(P,a  -P”)-V„i0a)(3Ptc  -P,A  -2P°)} 


+  i?r  In 


V vCH30H(l)A-H20(l)  J 


^  3/2\ 

p/,C 

c  -* 

no 


.  f5 

Tco,  r 


The  pure  component  molar  volumes  VCH  0H(1)  and  VIUO(l) 


(28) 

can  be  detennined  using  the 
methods  given  in  Thermodynamics  by  Smith  et  al.  (pp.  94-96). 

The  total  pressures  in  the  anode  and  cathode  plenums  may  be  expressed  as: 

pt,A  =  p4  +  pA  +  pA  (29a) 

1  rC02(g)  ^  rCH,OH(g)  ^  rH20(g)’  V-yA) 

the  partial  pressures  of  CH3OH,  Pch2oh(S)  ■>  H20, 
following  equation: 


Pi(g)  y«S)P  Xi(l)Pi 


^  0(  >  may  be  obtained  from  the 

(29b) 


(Raoult’s  law  for  the  ideal  solution  and  ideal  vapor  mixture  behavior) 


where  y;(g)  =  mole  fraction  of  species  i  in  the  gas  phase, 

Xi(i)  =  mole  fraction  of  species  i  in  the  liquid  phase, 

P°at  =  saturated  vapor  pressure  of  pure  i  at  the  cell  temperature  can  be  estimated  from  an 
empirical  equation  such  as  the  Antoine  equation. 


From  Eq.  (29a), 


1  =  Tco, 


<g) 


+ 


f  pA 

rCH,OH(g) 

r>t,A 


+ 


f  pA 
rH2Q{g) 
r>t,A 


(29c) 

Calculate  Vco,(g)  from  Eq-  (29c)  and  use  in  Eq.  (28). 

Also,  Pt  C  =  Pn20(g)  +  Po-Ag)  (f°r  the  case  of  pure  oxygen  as  the  oxidant)  (29d) 

or 

=  fn’oig)  +  po2(g)  +  pn2  (g)  (f°r  the  case  of  air  as  the  oxidant)  (29e) 

From  Eq.  (29d), 
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1 


f  pC 
1  H  2  0(g) 
D*,C 


_l_  c 

+yo2(S) 


(29f) 


From  Eq.  (29e), 


1  = 


H20(g) 

T)t,C 


+ 


f  Pc 

ro2{g) 


-,t,C 


J  V 


+ 


JV2(g) 

nt.C 


(30) 


In  the  gas  phase, 


Therefore, 


p(C) 

Uv2(g) 

p(C ) 

o2(g) 


P(C)  i 
Uv2(g) 


3.76 

1 


(for  the  large  excess  air  condition)  (31) 

(air  composition) 


3.7677 


(Ms) 


(32) 


Therefore, 


p(C) 

Uv2(g) 

T)t,C 


3  16P(C) 
J-,uro2(g) 


(33) 


Using  the  result  from  Eq.  (33)  in  Eq.  (30),  the  following  results: 


1  =  li^  +  4>76 


(  pC  \ 
ro2(g) 

pl.C 

V  ,7 


1  = 


Pc 

^O(g) 

D'.C 


+  4.76_y, 


02(g) 


Hence,  for  the  air  feed  to  the  cathode: 


(34) 
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or 


Also, 


1  = 


pc 

H2Q(g) 

nt,C 


c  c 

+  y02(g)  +  ->A2(g) 


(35) 


•  The  evaluation  of  the  standard-state  chemical  potentials  //(A’S  of  the  species  at  any 

temperature  T  other  than  their  values  given  at  a  reference  temperature  such  as  25°C, 
i.e.,  T0  =  298. 15K,  can  be  carried  out  by  the  thermodynamic  equations  developed  below: 

The  change  in  chemical  potential  of  a  pure  species  is  given  by 


d/ui  =  Vi  dP-  S;dT  (36) 

[Ref.  Smith  et  al.  Introduction  to  Chemical  Engineering 
Thermodynamics,  p.  181  (1996  edition).] 


If  the  standard  state  of  a  species  i  is  chosen  at  the  pressure  of  P  =  P°  =  lbar;  then, 

Eq.  (36)  is  written  for  a  selected  standard  state  of  a  species  at  the  constant  pressure  of  P  =  P° 
as 

d/uf  =  -SfdT  (37) 

The  change  in  the  chemical  potential  /uf  at  the  reference  pressure  P°  for  the  change  in 

temperature  from  T0  (reference  temperature)  to  any  temperature  T  is  obtained  by 
the  integration  of  Eq.  (37): 

ftA=/',A  T 

J  d/uf  =  -J  SfdT  (38) 


T 


T0 


▼ 

T 

=  AG  AfJo  -  J  SfdT  (at  P°=l  bar) 

Ta 


(39) 

(40) 

(41) 
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For  a  pure  chemical  species  i,  change  in  entropy  per  mole  for  a  change  in  its 
thennodynamic  state  is  given  by: 

,dT  (SV^ 


dS,  =  C„T- - 


\dT  jP 


dP 


(42) 


(p.  184  in  Ref.  Smith  et  al.  Thermodynamics  (1996 

edition)) 

where  Cp;  =  heat  capacity  per  mole  of  species  i, 

Vi  =  molar  volume  of  species  i  (@T,P). 

For  a  selected  standard  state  of  a  pure  species  i  at  the  reference  pressure  of  P°  (=1  bar),  Eq. 
(42)  is  reduced  to 


dS?=Ci* 

l  pi  rji 


(43a) 


Integrating  Eq.  (43a), 


S,  =  S*  + 


lc 


dT 

T 


(43b) 


where  SfT  =  entropy  per  mole  of  pure  (standard  state)  species  i  at  the  reference 


C 


pressure  P°  (=1  bar)  and  the  reference  temperature  To  (=  298. 15K  or  25 
°C), 

=  heat  capacity  of  pure  species  i. 


Insert  Sf  from  Eq.  (43b)  into  Eq.  (41)  to  obtain: 


-  \ 


dT 


Aa=AGV.-JV/7'-j 


1 

1 


rA  — 

^  pi  rp 


dT 


M, 


P^Gt,,T.-KTST-T,)-\ 


1 

1 


r&  dT 

^ PI  rj. 


dT 


(44) 
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where  A Gff  T  =  Gibbs  free  energy  of  formation  of  species  i  in  the  selected  standard 

state  (from  its  elemental  species  in  the  pure  state  at  the  standard  state 
pressure  of  P°=  1  bar,  and  the  reference  temperature  of  T  (=298.15  K), 
SfT  =  entropy  per  mole  of  species  i  in  the  standard  state  at  the  reference 
temperature  of  To  (=298.15  K), 

Cpj  =  constant  pressure  heat  capacity  of  pure  species  i  in  the  selected  standard 
state. 


The  value  of  SfT  is  to  be  computed  from  the  following  equation: 

.  A HffT  -  AG*  _ 

Kt.  =  \  +2>A7 


(45) 


where  data  on  A Hff  T  (enthalpy  of  formation)  and  AG; A  T  (Gibbs  free  energy  of  fonnation)  at 


the  reference  temperature  of  To  =  298.15  K  are  usually  provided  in  the  thermodynamic 
literature(e.g.  Smith  et  ah,  Introduction  to  Chemical  Engineering  Thermodynamics  (1996 
edition)),  o  j  is  the  number  of  moles  of  an  elemental  species  j  involved  in  fonning  1  g-mole  of  a 


chemical  compound  i  and  Sf  is  the  standard  state  entropy  per  mole  of  the  elemental  species,  j, 


at  the  reference  temperature  of  T0  (=298.15  K).  Data  on  the  standard  state  entropies  of  elemental 
species,  S'f  ,  are  also  provided  in  the  literature  (e.g.  R.C.  West  (editor),  Handbook  of 


Chemistry  and  Physics,  CRC  Press,  Cleveland,  Ohio.) 

The  standard  state  heat  capacity  of  a  species  i  in  the  gas  phase  at  temperature  T  is  given 


as:  Cf=R 


A  +  BT  +  CT2  +^r 


(46) 


where  R  =  8.3 14  J  mol"1  K"1  (the  universal  gas  constant) 

T  =  temperature  in  [K] 

A,B,C,D  =  heat  capacity  coefficients  for  the  species  i  in  the  ideal  gas  standard  state  (values 
given  on  p.  638  in  the  reference  given  above). 


The  standard  heat  capacity  of  species  i  in  the  pure  liquid  phase  at  temperature  T  is  given 
by;  C%  =r[a  +  BT  +  CT2]  (47) 

where  values  of  the  heat  capacity  coefficients  A,B,  and  C  for  the  pure  chemical  species  in  the 
liquid  phase  are  given  on  p.  639  in  the  reference  by  Smith  et  al. 

The  molar  volume  of  a  chemical  species  i  in  the  liquid  phase  at  the  mixture  liquid  phase 
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temperature  and  pressure  may  be  computed  from  the  following  scheme: 

The  Rackett  equation  can  be  used  to  obtain  the  molar  volume  of  a  pure  liquid  in  the  saturated 
state,  V- ,  at  the  temperature  T  of  the  fuel  cell  system: 


V  sat  =  V  7 

v  i(l)  yci^c 


(l -Trf 


where 


IW 

T„(K) 


(reduced  temperature) 


(48) 


Vci ,  Zci ,  Tci  =  critical  molar  volume,  compressibility  factor,  and  temperature  of  the 

species  i.  The  species  critical  property  data  are  given  for  some  chemical  species  in  Appendix  B 
in  the  reference  by  Smith  et  al  (1996  edition  [pp.  634-636]).  Also,  one  detennines  the  saturated 
vapor  pressure  of  the  pure  liquid  phase  species  i  at  the  system  temperature  T  using  the  following 


p 

Antoine  equation:  log10  P;sm  =  A-  — — — 

where  P-at  =  vapor  pressure  of  pure  i  in  the  liquid  phase  at  the  temperature  T  (in  °C) 


(49) 


A,B,C  =  the  Antoine  equation  constants;  values  of  which  are  given  in  the  Ref.  R.M. 

Felder  and  R.W.  Rousseau,  Elementary  Principles  of  Chemical  Processes,  pp. 
640-641,  John  Wiley  &  Sons,  Inc.  (2000,  third  edition). 

psat 

Then,  one  determines  the  reduced  density,  ps‘n  at  P/f  =  — —  using  the  generalized 


density  correlation  for  liquids  given  on  page  95  in  the  Ref.  by  Smith  et  al.  Also,  at  the  system 

p  T 

pressure  and  temperature,  using  Pr=  —  and  Tr  = — ,  one  determines  pr.  Then,  the  required 

ci  ci 

molar  volume  Vi(i)  of  a  pure  species  in  the  liquid  phase  is  computed  from  the  following  equation: 

-p  sat 

Vm(@T,P)  =  psrat^  (50) 

Pr 

Comments: 

•  The  Rackett  equation  is  accurate  within  5%  [See  the  Ref.  by  Smith  et  al.]. 

•  Because  the  effect  of  pressure  on  the  liquid  phase  molar  volume  is  small  for  pressures 

<  15  bar;  then,  one  may  assume  ViU)[@T,P]  -  ( @T )  for  pure  species  i. 
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4.1. 1.2  Calculation  Scheme 

A.  For  CH30H(1)  in  FEC^l)  feed  composition  at  a  temperature  T,  calculate  the 
reversible  or  open-circuit  cell  voltage,  U  for 

PtC  =1,3, 5, 7  bar;’ 

for  Pt,c  and  Pt  A 

Pl’A  =  1, 3, 5, 7  ba 

to  be  of  equal  and  different  values.  Repeat  such  calculations  for  temperatures  of  25,  50,  60  and 
80°C  for  each  fuel  feed  composition. 

B.  For  the  case  of  fuel  feed  composition  of  ^  1  molar  CH3OH(l)  in  water(l),  one  may 

use  Eq.  28  and  for  the  case  of  fuel  composition  of  >  1  molar  CH30(1)  in  water(l),  one  should  use 
Eq.  (27).  In  case  the  total  pressures  in  the  anode  and  cathode  plenums  are  less  than  10  bar,  one 

may  set  =  1,  i.e.,  the  gas  phase  is  assumed  to  follow  ideal  behavior. 


4.1. 1.3  Sample  Calculations 


[Note:  This  computation  is  only  for  the  illustration  of  the  calculation  steps.] 

(1) 

•  Cell  temperature,  T  =  25°C  =  298. 15K 

Fuel  feed:  1  molar  CH3OH  in  FEO(l)  at  Pt  A  =  1  bar. 

•  Air  feed:  Pt  C  =  1  bar. 

P°  =  1  bar  -  standard- state  pressure. 

From  the  thermodynamic  tables,  @T  =  25°C  =  298. 15K: 


A  _  A 

Ach ,OH(l)  ~ 


McH,OH(l) 


■  166270 J 
mol 


’  Ao 2(g) 


=  0.0 


J 

mol 


Aco2(g)  =  -394359- 


Ah2o(1)  =  -237129- 


R=8.314  J/mol-K 
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1  g  —  mol 

In  one  molar  methanol  solution  of  the  anode  feed,  nCH  0H(l)  =  - ;  and 


1  liter* 


gm 

\ 

_  I  _ 

) 

1  g  -  mol  *  M 


CH,OH 


l 


gm 

g  -  mol 


1  H20(l) 


M 


(51) 


h2o 


(where  psolu  =  density  of  the  solution) 


vCH3OH(l) 


n 


CH3OH(l) 


nCH,OH(l)  +  nH20{l) 


A  _  I  _  A 

XH2OU)  —  1  —  XCH3OH(l) 


(51a) 

(51b) 


P water  (I) pure  25  C- 


1  kg 

\m3 

=  0.997 

\kg] 

_  0.00 1003m 3  _ 

1000/ 

_  1  _ 

(frorn  steam  tables  p.  662-Felder  and  Rousseau) 


nna(l)  ( estimated )  = 


( 

0.997  — 

(i/) 

^lOOOgm^ 

lg  -  mol  ^ 

LI  i  J 

l  lkg  J 

L  18  gm  )_ 

=  55  Ag  -  moles 


(See  Sections  C  and  D  for  the  exact  and  approximate  procedures  to  compute  the  liquid  phase 
mole  fractions,  pp.  40  -  49) 

1  g  -  mol 


vCH3OH 


(l  +  55.4)g  -mol 


=  0.017731 


xh2oo)  —  1  xch,oh(i )  ~  0.982269 


-VC02(g )  ? 


(51c) 
(5  Id) 


From  Eq.  (29c) 


l-yCo2(g)  + 


f  pA 
rCH3OH(g) 


a  r  pa  \ 

rH20(g) 


)t,A 


V 


+ 


■)t,A 


From  Eq.  (29b), 

A 


pA  _  A  psat 

rCH3OH(g)  aCH3OH(1VCH3OH 


pA  _  A  psat 

rH20(g)  ~XH20(l)rH20 


at  25°C  from  the  Antoine 
Eq.  (49) 


(52) 


(52a) 


(52b) 
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1473  11 

lo§i0  Pch-soh  =  7.87863  -  2  (vapor  pressure  in  mm  Hg) 

=  2.101728 

(A,B,C  data  from  Felder  and  Rosseau) 


log10P^  =8.10765- 


1750.286 
25  +  235 


(vapor  pressure  in  mm  Hg) 


=  1.375781 


Pch3oh  =  126.4mm Hg  ;  P^0  =  23.76mm  Hg 


Pchoh  =  (126.4 mmHgi  LQ1325^r 
C"P"  V  *  1  760mm  IIg  j 


=  0.1685bar 


PZ‘0  =  (23.76mm  Hg) 


A  1.0 1325 ba/P 


=  0.0317  bar 


y  760mm  Hg  j 

All  the  parametric  infonnation  into  Eqs.  (52a)  and  (52b): 
Pc„}OH(g)  =  0.017731*0.1685  bar  =  0.00299  bar 


H20(g) 


=  0.982269  *  (0.03 17  bar )  =  0.03 1 14  bar 


Substitute  the  data  above  into  Eq.  (52)  with  P,  A  =  P"  =  1  bar 
1  =  ycoAg)  +  (0.00299)  +  (0.031 14) 


Aw  =0-9659 


Now>  y  o2(g) = 2 
From  Eq.  (34a), 


c 

yo2(g) 


4.76 


1- 


pc 

mo(g) 

D  t.C 


Using  Eq.  (29b): 


pL  _  C  psat 

rH20{g )  ~  XH20{l)rH20 


25°C 


PH2o(g)  =0.0317  bar 
From  (52f),  (52h),  with  Pt-C  =  lbar: 


c 

-vo2<g) 


^j[l- 0.03 17]  =  0.203 


(52c) 


(52e) 


(52f) 


(52g) 

(52h) 


(52i) 
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Insert  all  the  parametric  data  into  Eq.  (28)  to  obtain: 


Term  I 


FU  = 


-166270  + 1  (0)  -  (-  394359)  -  2(-  237129); 


J 


[(  Vch3oh(0  0  - 1,  bar)-  VHi0{f)  (3  *  1  - 1  -  2, bar)]1 


mol  -  CH3OH  j 


8.3147 
mol  -  K 


(298.15^)ln 


3 


(.01773 1)(.982269) 


0.203 


1  bar 
Ibar 


0.9659 


1  bar 
Ibar 


JJ 


i 

L. 

Term  III 

Term  II 


FU  =  f  ^  j[(702347)  +  0  +  (- 1 5,882.934)] 


Term  I 

Term  II 

Term  III 

J 

mol 


F  =  96487^?^  =  96487-  C°“hmb 


U- i 
6 


equiv 

702347 


96487 


■  + 


0 


96487 


+ 


1  mol  of  unity  charged  particles 
J 

mol 


-15882.934 


V 


96487 


J 


C 

mol 


Note:  [J  =  Coulomb'Volt] 
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U  =  -  [7.2792  +  0  +  (-  .1646)]  Volt 


•  Cell  open-circuit  potential  for 

•  Pt,A  =  Pt  C=l  bar,  air  as  oxidant, 

•Fuel  composition:  1  molar  CH3  OH  (in  water  (1)), 
•Cell  temperature  T=25°C 


(53) 


Term(llf) 

-0.023 

U 

1 

2.3% 


Then,  @  T=60°C;  P 


t,A 


'  <3% 

=  Pt  C  =  lbar;  1  molar  feed  of  CH3  OH(l)  in  H20(1),  air:  oxidant 


L 


First  calculate  juf's  at  T=60°C; 
Then,  calculate  U 


Here,  activity  coefficient  values  would  be  needed!  (one  should  make  calculations) 

4.1. 1.4  Activity  and  Fugacity  Coefficients 

The  local-composition  models,  though  limited  in  flexibility  in  the  fitting  of  experimental  data, 
are  adequate  for  most  engineering  purposes.  These  can  be  generalized  to  multicomponent 
systems  and  require  only  the  parameters  needed  to  describe  the  constituent  binary  systems. 
Here,  the  Wilson  equation  for  a  multicomponent,  liquid  phase  has  been  used  and  is  given  by 


— =-Z 

RT  “ 


( 

\ 

x,  In 

j_ 

(54) 
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(55) 


l  J  J  *  E*A, 

j 

where  GE  =  G  -  Gldeal  solutlon  =  excess  Gibbs  free  energy  per  g-mole  of  the  liquid  phase  at  the 
liquid  phase  at  the  liquid  phase  temperature,  pressure,  and  composition  conditions;  x,  (or  xj), 
mole  fraction  of  a  component  in  the  liquid  phase;  y.  ,  the  activity  coefficient  of  a  chemical 
species  i  in  the  liquid  phase.  In  the  above  equations  all  summations  are  over  all  chemical 
species.  For  each  ij  pair,  there  are  two  interaction  parameters  such  as  A.,  and  A  .. ,  because 

A..  ^  A  /V .  Fori=j,  A(y  =  l.  The  temperature  dependence  of  the  parameters,  A..,  is  computed 
from  the  following  equation: 

Aij=17  expf *  j) ’  (56) 

vi  V  R1  J 

where  Vj  and  V;  are  the  molar  volumes  at  temperature  T  of  pure  liquids  j  and  i,  and  a,,  is  a 
constant  independent  of  composition  and  temperature.  Values  of  V;  and  ay  of  some  binary 
systems,  including  CH3OH  /  water,  are  given  in  Smith  et  ah,  Introduction  to  Chemical 
Engineering  Thermodynamics,  p.  467,  The  McGraw-Hill  Co.  (1996). 

Fugacity  coefficient  of  a  component  k  in  a  multicomponent  gas  phase  is  computed 
from  the  following  general  equation: 


the  dummy  indices  i  and  j  run  over  all  species  in  the  gas  phase  mixture,  and 

4=2 (58a) 

4=2  4-4-4’  (58b) 

with  8U  =  8kk  =  0 ,  etc.,  and  8 ki  =  8ik ,  etc. 

Values  of  the  viral  coefficients  of  pure  chemical  species  B;;,  By,  etc.,  can  be 
calculated  from  the  following  generalized  correlations: 


=  B°  +  coB' 


(59) 


The  second  viral  coefficients  are  functions  of  temperature  only.  Similarly  B°  and  B1  are 
functions  of  reduced  temperature,  Tr  only. 
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B°  =0.083 


0.422 


(60) 


t1.6 


B 1  =  0.139  — 


T 

r 

0.172 

rp  4.2 


(61) 


The  cross  coefficients  B,k,  By,  etc.,  are  computed  from  the  following  equations: 

(62) 


RT 

*cij 


where  B°  and  B 1  arc  the  same  functions  as  given  in  Eqs.  (60)  and  (61).  The  combining  rules  are 
used  to  compute  colt ,  77. ,  Pcj  as  follows: 


CO;  +6): 

coji  = - - 

■J  2 


(63a) 


=  tM-Kj) 


z  +z 

ci  cj 


v«j  = 


vT+vf^ 


PdJ  = 


ZcyRTcli 

z. 


(63b) 

(63c) 

(63d) 

(63e) 


In  Eq.  (63b),  ky  is  an  empirical  interaction  parameter  associated  with  an  i-j  molecular  pair. 
When  i=j  (i.e.  a  pure  chemical  species)  and  for  species  that  are  chemically  similar  ky  =  0.  In  the 
absence  of  data,  ky  is  set  =  0. 


The  reduced  temperature,  Tnj 


T 

- ,  for  a  given  i-j  pair,  is  used  for  the  computation 

Pcij 


of  B°  and  B1  from  Eqs.  (60)  and  (61)  for  application  in  Eq.  (62).  For  more  information  on  the 
subject,  Molecular  Thermodynamics  of  Fluid-Phase  Equilibria  by  Prausnitz  et  al.  (1999)  may  be 
consulted. 


In  the  above  equations  coi  and  77 ,  Zci ,  and  Vci  are  the  accentric  factor  and  critical 

temperature,  compressibility  factor  and  volume  of  a  species  i  respectively.  Values  of  these 
parameters  are  given  in  the  Ref.  by  Smith  et  al.  (1996),  pp.  635,  636. 
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4.1. 1.5  Discussion  of  the  Computed  Results 

Set  1  -  Methanol  Concentration  Range:  0.01-0.50  M  (Figs.  2  through  15). 

Figures  2  and  3,  and  9  and  10  show  the  effect  of  pressure  on  the  cell  voltage  at  the 
cell  temperature  of  25  °C  at  three  molar  concentrations  of  methanol  in  water  as  the  fuel  for  the 
ideal  and  nonideal  liquid-phase  solution  and  gas  behavior,  respectively.  It  is  quite  apparent  that 
for  each  fuel  concentration  the  open-circuit  voltage  increases  with  an  increase  in  the  cathode 
pressure  over  the  range  of  1-10  bar,  whereas  it  decreases  with  an  increase  in  the  anode  pressure 
over  the  same  pressure  range.  Figures  2  and  9  show  the  effect  of  accounting  for  the  nonideal 
behavior  of  liquid  and  gas  phase  on  the  open-circuit  voltage  as  function  of  cathode  pressure  and 
fuel  concentration.  There  is  a  slight  increase  in  open-circuit  voltage  in  the  case  of  accounting  for 
the  nonidealities  of  liquid  and  gas  phases.  Similar  effect  of  the  nonideal  behavior  of  liquid  and 
gas  phases  is  shown  by  Figs.  3  and  10  at  each  anodic  pressure  and  fuel  molar  concentration. 

Figures  4  through  8  show  the  effect  of  temperature  on  the  open-circuit  voltage  for  the 
ideal  liquid-phase  solution  and  gas  phase  behavior.  The  voltage  decreases  with  an  increase  in  the 
cell  temperature  at  a  set  of  fuel  concentration,  cathodic,  and  anodic  pressure  values.  Similar 
trend  is  seen  in  Figs.  11  through  15  for  the  nonideal  liquid-phase  solution  and  gas  phase 
behavior. 

The  overall  conclusion  from  the  computed  data  is  that  a  set  of  higher  cathode  and 
lower  anode  pressures,  higher  methanol  concentrations  in  the  fuel,  and  lower  fuel  cell 
temperature  results  in  a  higher  open-circuit  voltage.  For  example,  at  the  cathode  pressure  of  10 
bar,  anode  pressure  of  1  bar,  methanol  concentration  of  0.5  M  in  the  liquid  fuel  and  temperature 
of  25  °C  with  air  as  the  oxidant  in  the  cathodic  chamber;  the  open-circuit  voltage  is  more  than 
1.200  volt  shown  in  Fig.  13;  in  comparison  with  less  than  1.145  volt  at  the  cathode  pressure  of  1 
bar,  anode  pressure  of  10  bar,  0.01  M  methanol  concentration,  and  temperature  of  60  °C  as 
shown  in  Fig.  8. 

Figures  16  through  29  show  the  computed  results  for  the  assumption  of  water 
formation  at  cathode  in  vapor  phase  for  the  Set  1  -  Methanol  concentration  range:  0.01  -  0.50  M. 
Some  effect  of  this  assumption  is  shown  by  the  data  in  Figs.  25  through  29  relative  to  the  data  in 
Figs.  1 1  through  15  for  the  situation  of  water  fonnation  at  cathode  in  the  liquid  phase.  However, 
one  observes  the  similarity  in  the  trends  of  the  plotted  data  for  these  two  situations. 
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Sets  2  and  3  computed  data  are  for  methanol  concentration  ranges  of  0.5  -  1.0  and  1.0 


-  2.0  M.  Trends  of  the  plotted  data  for  these  sets  are  similar  to  Set  1  data. 

4.1.2  Equilibrium  Open-Circuit  Voltage  Equation  Adjusted  for  the  Case  of  Water 
Formation  in  the  Vapor  Phase  by  Hydrogen  Ion  Oxidation  at  the  Cathode 

(see  sketch  on  page  5) 

The  cell  reactions  are  written  as: 

At  anode: 


CH3OH  (1)  +  H20  (1)  ^  C02  (g)  +  6H+  (1)  +6e 


(64) 


( in  |3  phase) 


At  cathode: 


3/2  02  (g)  +  6H+  (1)  +6e  ^  3H20  (1) 


3H?Q  (1)  <->  3Tf>0  (g) 


3/2  02  (g)  +  6H+  (1)  +6e  <->  3H20  (g) 


(65) 


The  overall  cell  reaction  is 


CH3OH  (1)  +  H20  (1)  +  3/2  02  (g)  <->  C02  (g)  +  3H20  (g) 


(66) 


The  cell  potential  equation  is 


6 


f 


+  RTlnt 


V  CH2OH(l)Y 


H20(l)XCH20H(l)XH20(l) )  , 


W 


iA  A  r 

rco2yco2  p, 


(67) 
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For  the  ideal  liquid  phase  and  gas  phase  behavior, 


FU 


l 

6 


A  A  .  ^  A  A  'j 

" CH,OH(l )  "c02  +  T/‘o2  +  "h20(D  ~  ^ M 


H20(g) 


(VcH,OH(n  +  Vh20(1)\p  '  P  ) 


J 


+  RTln 


\XCH30H(I)XH20(l)  Ji 


(  pt,C  \ 

c  £_ _ 

y  o2  po 

V  r  J 


3/2 


It, A  \ 


A  1 

yco2  o 
W  r  J 


*,C  \ 


c 

y h,o 


J  J  J 


(68) 


Standard-state  potential  U9  from  Eq.  (68): 


U 9  = 


A  _  A  i_,,A  .  A  _ -j  A 

MCH3OH(l)  Mc02(g)  +  ^  Ao2(g)  +  Mh20(1)  ^^H20(g) 

6F 


(69) 


Using  the  thermodynamic  equations  given  in  the  document  A.,  Eq.  (69)  has  been  reduced  to 


U 9  = 


A  _  A  A  _i_A  _  o  A 

McH,OH(l),T0  Mc02(g),T0  +  ^Mo2(g),T0  +  Mh20(J),T0  ->/UH20(g),T 

6F 


+{t-to\ 


T. 

_  oA  ,  oA  _ _  cA  _  cA  ilOA 

°CH2OH0),To  °C02(g),To  2  J H20(l),T0  "r  H20(g),T0 

6F 


(70) 


( =  ,  approximate  because  heat  capacity  correction  terms  have  been  neglected).  T0  =  298. 15  K, 
the  reference  state  temperature. 


4.1. 2.1  Sample  Calculations 

Using  the  thermodynamic  data,  the  following  approximate  equations  for  U9  have  been  obtained: 
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U9  =1.2132-  0.00 14(T-T0),  [vo/t]  (water  produced  at  the  cathode  in  liquid  phase)  (71a) 
U9  =1.1 689  +  0.00048(7"  —  To  ),  [vo//]  (water  produced  at  the  cathode  in  gas  phase)  (71b) 


Temperature  [K] 

Ue  from  Eq.  (71a) 
[voltl 

Ue  from  Eq.  (71b) 
[voltl 

323.15 

1.1781 

1.1809 

353.15 

1.1362 

1.1953 

Comment:  Note  the  effect  of  temperature  on  approx.  Ue  .  If  the  water  produced  at  the  cathode 


is  taken  to  be  in  the  liquid  phase,  Ue  decreases  with  an  increase  in  temperature.  However,  it 
increases  with  an  increase  in  temperature  if  the  water  produced  at  the  cathode  is  taken  to  be  in 
the  gas  phase. 

4.1. 2.2  Actual  Procedure  to  Determine  the  Liquid  Phase  Mole  Fractions  of  CH3OH 
and  H20,  given  the  Molarity  of  the  Solution 


In  general,  the  liquid  phase  mixture  containing  ncil0ll  moles  of  methanol  and  nHO  moles 


of  water  at  given  pressure  and  temperature  values  has  the  volume  given  as: 


n  v  -| -n  V  -V' 

nCH,OH  r  CH2OH  T  nH20  r  H20  r  solution 


(72) 


where  VCH  OH  and  V n  0  are  the  partial  molar  volumes  of  methanol  and  water  in  the  liquid 


mixture;  expressed  as  liter  per  g-mol  of  a  species.  By  definition,  methanol  molarity  is, 


CM  = 


,'CH2OH 


V. 


solution 


[mol  of  CH3OH  per  liter  of  solution] 


(73a) 


Also, 


lH-,0 


V1 

solution 


(73b) 


Mole  fractions  of  methanol  and  water  are  then  given  as 


lCH,OH 


'  M 


vCH,OH 


n  CH-.OH  +  nH20  CM  +  CW 


XH,0  — 


1  HnO 


nCH,OH  nH.O 


CM  +  CW 


(73c) 

(73d) 
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where  xCH^0H  +  xHi0  =  1.0 


(73  e) 


Equation  is  written  as 


,  v‘ 

V  solution  J 


V  + 

v  CH-.OH  ^ 


,1h2o 

.  V 

V  solution  J 


VHl0  =  1 


c  V  +r  V  =1 

CH.OH  ~  ^Wv  H.O  1 


(74) 


(74a) 


If  the  molarity  of  the  solution  with  respect  to  methanol  has  been  fixed,  then  one  should 
determine  cw  ,  molar  concentration  of  water  in  the  liquid  mixture  to  calculate  the  mole  fraction 


of  methanol,  xCH^OH  ,  from  Eq.  (73c).  This  can  be  done  with  the  help  of  Eq.  (74a)  if  one  has  the 

values  of  partial  molar  volumes  of  the  species  at  the  desired  temperature  and  pressure  conditions. 
The  partial  molar  volumes  are  functions  of  the  species  mole  fractions  at  the  given  temperature 
and  pressure  conditions.  These  should  either  be  determined  experimentally  (using  the  definition 
of  a  species  partial  molar  volume)  or  from  the  literature  data  -  in  the  graphical  or  mathematical 
equation  form.  If  one  likes  to  use  the  data  on  the  species  (CH3OH,  H2O)  partial  molar  volumes 
as  given  in  the  graphical  form  in  Fig.  10.3  of  the  reference:  J.M.  Smith,  H.C.  Van  Ness,  M.M. 
Abbott,  Introduction  to  Chemical  Engineering  Thermodynamics,  p.  327,  the  McGraw-Hill 
Companies,  Inc.  (1996,  5th  edition);  one  should  use  the  following  procedure: 

First  assume  that  the  mixture  follows  the  ideal  solution  behavior,  then  VCH  OH  =  VCH  OH , 

VH  0  =  VH  0 ;  i.e.  partial  molar  volumes  are  equal  to  the  pure  component  molar  volumes  at  the 
temperature  and  pressure  conditions  of  the  desired  solution.  Equation  (74a)  is  written  as: 
CM^CH,OH  +CW^H20  (74b) 

For  a  fixed  cM  value,  values  of  pure  component  volumes  Vcn  o,,  and  VH  O  at  the  solution 
temperature  and  pressure  are  inserted  in  Eq.  (74b)  to  obtain  the  value  of  cw  .  The  mole  fraction 


of  methanol,  xCH  On  ,  should  then  be  computed  from  Eq.  (73c).  At  this  composition  the  partial 


molar  volumes  are  obtained  from  Fig  10.3  and  cw  is  now  computed  from  Eq.  (74a)  for  the  fixed 
cM  value.  Calculation  of  xclf  OII  is  repeated.  The  most  recently  calculated  value  ofxc//  0//  is 
compared  with  its  previous  value.  If  these  values  agree  within  a  preset  tolerance;  then,  the  most 
recently  calculated  values  of  cw  and  xCHOH  for  a  fixed  cM  value  at  the  desired  temperature  and 


35 


pressure  conditions  are  correct.  If  the  most  recent  and  previous  value  ofxCH  OH  do  not  agree 


within  the  set  tolerance,  the  procedure  is  repeated  until  the  tolerance  condition  is  satisfied. 
Example  1 .  What  is  the  mole  fraction  of  methanol  in  the  1  molar  solution  of  methanol  in  water 
at  25  °C  and  1  atm. 

Solution: 

cM  =  \mol/ liter  of  methanol  in  the  methanol-water  liquid  at  25  °C,  1  atm. 

At  25°,  1  atm,  the  pure  component  molar  volumes  are 
VCHj0h  =  0.040727  liter/ mol ;  Vu_0  =0.018068 liter / mol . 


From  Eq.  (74b): 

_  1 —  cm  ^ch3oh  1  -  (l)(0. 040727)  _  mol  of  water 

Cw~  Vf0  ~  0.018068  ~  '  liter 


'/CH,OH 


CM  CW 


- =  0.01849 

1  +  53.092 


(75) 

(76) 


From  Fig.  10.3  in  the  reference  by  Smith  et  ah,  at  xCHOH  =  0.018149; 


VCIIOH  =  V<zh,oh  =0.037 5  liter /mole  and  VHi0  =  VH  O  =0.01843  liter /mole 
From  Eq.  (74a): 


1-c  V 

1  r  CH-,OH 


V , 


1  -  (l)(0.0375)  _  23  271  mo^  °f  water 


H.O 


1 


lCH3OH 


1  +  53.271 


0.018068 


=  0.01843 


liter 


This  value  differs  from  that  in  Eq.  (76)  by 
10.01843  -  0.018491 

- - -x  100  =  0.326%  [i.e.  <  1%  (set  tolerance)] 

0.01843 


(77) 

(78) 


Repeat  the  procedure. 

At  xcll  0ll  =0.01843,  VCHiOH  and  Vu0  are  obtained  from  Fig.  10.3.  They  are  about  the 


same  as  given  above.  So,  xCH^OH  =  0.01843  as  in  Eq.  (78).  Thus,  no  further  trial  is  needed. 


xCHiOH  =0.01843  =  0.018,  xHi0  =1-0.01843  =  0.98157  ^0.982 

On  the  question  of  effect  of  pressure  on  the  species  partial  molar  volumes  in  the  liquid 
phase  mixture  of  methanol  and  water,  one  may  ignore  the  effect  of  pressure  on  the  molar 
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volumes  over  the  pressure  range:  1-10  bar.  This  is  supported  by  the  information  given  on  page 
185  of  the  reference  by  Smith  et  al.  For  example,  molar  volumes  of  pure  liquid  water  at  25  °C 
are  18.071  and  18.070  cm  /mol  at  pressure  of  1  and  10  bar,  respectively.  As  to  the  effect  of 
temperature  on  the  partial  molar  volume  of  a  species  (CH3OH,  H20)  in  the  liquid  mixture,  one 
may  obtain  the  parameter  /?. ,  from  the  experimental  data  on  the  partial  molar  volumes  as  a 


function  of  temperature  at  fixed  pressure  and  composition  conditions,  defined  below: 

'aFA  M  (79) 


a  4 


V  dT  jp 


where  Vl  =  partial  molar  volume  of  species  i  in  the  liquid  mixture  of  given 


composition,  temperature  and  pressure  conditions. 

Data  on  /?  versus  temperature  should  be  generated  as  function  of  temperature  for  each 


fixed  set  of  composition  and  pressure  conditions.  For  dilute  solutions  of  methanol  in  water  (i.e. 
for  cM  <  1  mol /liter  ),  /3W  for  water  may  be  approximated  by  the  pure  water  expansion 


coefficient,  (5W  =  — 


V, 


w  V 


dV 

'JV  w 
dT 


7  p 


The  experimentally  developed  infonnation  on  can  then  be  used  to  obtain  a  value  of  Vt 


at  the  desired  temperature  of  the  liquid  mixture  of  interest  provided  a  value  of  Vi  at  a  nearby 


temperature  is  known.  Note  that  the  partial  molar  volumes  of  the  species  are  required  in  Eq. 
(74a)  to  relate  the  species  molar  concentrations  and  mole  fractions  in  the  liquid  phase  mixture. 
Example  2:  Find  the  partial  molar  volume  of  water  at  50  °C  for  0.1  molar  solution  of  methanol 
in  water  at  pressure  of  1  bar. 

Solution: 

This  solution  is  very  dilute  with  regard  to  methanol.  At  25  °C, 

Vh2o,2s  =  Vh2o,  25  =  18.068  ctnmol1  (80) 


From  Eq.  (79), 


Pw  = 


1 


V, 


H.O 


dV, 


H,0 


1U  V 


dT 


J  P,x 

dVHi0  =(pwVHi0)dT 


(81) 

(82) 
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vH20,T  1 

J  av„t0  =  \(pwv„i0)dT 

yH2ojo  To 

T 

K,o.r  -yH,o,T.)=  l(firyH,o)dT 


l’ ll j  —  Vh,0.t„  4“  f  [ftw  Vh,o  'jdT  (at  given  composition  and  pressure  conditions)  (84) 

T0 

If  mathematical  expressions  for  fiw  and  VHO  are  available  as  functions  of  temperature, 
the  integration  can  be  carried  out  analytically  or  graphically.  If  such  information  is  not  available; 
then,  pw  and  VH  O  may  be  approximated  by  the  pure  water  expansivity,  (5W ,  and  pure  water 
liquid  volume,  VHO .  This  approximation  is  valid  particularly  in  the  dilute  solution  of  methanol 
in  water  of  0.1  molar  concentration.  For  such  a  situation,  Eq.  (84)  becomes 


T 

Vh20,T  ~  Vh20,T0  +  J"  ( ’ PwT H2o)dT 


One  may  use  the  arithmetic  average  of  the  ( PwVHi0 )  values  at  the  two  temperatures  in  the 
integration  on  the  right  side  of  Eq.  (85)  and  approximate  the  integral  as  given  below: 

Vh2o,t  =  Vh2o,t0  +  ( PwVH2o )(j  ~To)  (86) 


,  (o  T,  )  H20\iyc)  '  TW  H20\5()'CJ 

where  \pwVHi0)  = - - - 

_  [(256x  10  6 )(l 8.07 1)]+  [(456 x  10~6)(l 8.234)] 


+  {p> 


W  H2°  1 50-0 


,  \cm3  mol  lK  ‘J 


0.006470  I  cm  mot1  K  1 


Hence,  VH  OtSO  =  VH  Of2S  +(0.00647 cm3 mol  )([50  +  273. 1 5] -  [25  +  273. 15], 


18.068  +  (0.00647)(25)  =  18.230  cm3  mol 
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Note  that  this  volume  differs  from  the  given  value  of  VH  0  50  =  18.234  cm3 mol  1 
118.234 -18.2301 

by  - - -  x  100  =  0.022%  which  is  less  than  0. 1  %. 

18.234 

4.1. 2. 3  An  Approximate  Procedure  to  Determine  the  Liquid  Phase  Mole  Fractions 
of  CH3OH  and  H20,  given  the  Molarity  of  the  Solution 

If  the  solution  density  psolu  as  \gm  ■  cm  3  ]  is  available  as  a  function  of  the  molarity  of 
methanol  in  water,  cM  [g-mol  CH3OH  per  liter  of  the  solution]  at  a  given  set  of  temperature  and 
pressure  conditions;  then,  molar  concentration  of  water  in  the  solution  can  be  computed  from 


Psolu 

(  gms  ' 

f  1000 cm3'. 

- 

/  g  -  mol  Cl  1,0/1 ' 

(  gm  Cl fOfl  ) 

g  -  mol  H20 

l  3 

V  cm  J 

v  liter  y 

v  M  liter  y 

{  a,i°"  g-mol  Cl  If)  1 1 ) 

I  \ 

Mh2o 

gmH20 
g  -  mol  H20 

(90) 

Then,  xai/)ll  =  M 

CM  +  CW 

(90a) 

and  xHi0  =  1  —  xcll()ll 

(90b) 

In  the  event  that  neither  the  partial  molar  volumes  of  the  species  nor  the  solution  density 
as  a  function  of  composition  is  available;  then,  the  pure  component  methanol  molar  volume  at 
the  solution  temperature  and  pressure  is  used  to  approximate  the  volume  occupied  by  methanol 
in  the  solution  as  given  below: 

Volume  occupied  by  methanol  in  the  solution  of  cM  (mol / liter) ,  molar  concentration  of  methanol 
in  water,  is 


"  mol  CH3OH  ' 

^  V ch3oh’ l‘ter  °f  Pure  Cl I  .OI 1 3 

v  liter  of  solution  y 

v  mol  CII  fJH  y 

vm  =  CM  Vch,oh  \liter  of  CH2OH /liter  of  solution \  (91a) 

Therefore,  volume  occupied  by  water  in  solution  is  approximated  by: 

vwater  =  (l-vM)[liter  of  water / liter  of  solution ]  (91b) 

Now,  moles  of  water  present  in  1  liter  of  the  solution  are  approximated  by: 
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^  ^  liter  of  water 

liter  of  solution 

mol  H20 

liter  of  water 
mol  H20 

{  Vw  J 

liter  of  solution 

(91c) 


where  Vw  =  pure  component  water  molar  volume  at  the  solution  temperature  and 
pressure  conditions. 


Then,  the  mole  fraction  of  methanol  in  water  is  approximated  by: 


1 


'/CH,OH 


CM  +CW 


1  + 


\CM  J 


(9  Id) 


Using  Eqs.  (91a)  and  (91b), 


( c  ^ 

1 

1 

(  1  CM  VcH2OH  ^ 

= 

(  1  V  ^ 

1  v  CH^OH 

\CM  J 

CM 

V  Vw  ) 

CM 

l  vw  J 

CM  Vh20  VHio  j 

(91e) 


From  Eqs.  (6d)  and  (6e): 


x 


CHJOH 


r  1  V.  A 


1  + 


CH,OH 


\  CM  Vh20  Vh20  J 


CM  Vh^O 


CM  Vh,0 


c  V  + 1  -r  V 

r  H^O  T  1  r  CHJOH 


1  +  CM  Yh20  VcH,OH  ) 


(9  If) 


where  cM  = 


mol  of  Cl I 3  OH 

V  - 

liter  of  pure  CII  f)I I 

liter  of  solution 

5  v  ch3oh 

g-mol  of  CH3OH 

V, 


h2o 


liter  of  pure  H20 
g-mol  of  H20 


Equation  (9 If)  is  the  approximate  formula  to  calculate  xCH  OH . 

Comparison  of  the  actual  procedure  (Section  C)  with  the  approximate  Eq.  (9  If)  to  compute 
xCH  OH  given  the  molar  concentration  of  methanol  in  water: 


c 


M 


1  mol  of  Cl  If )I  I 
1  liter  of  solution 


at  25  °C,  1  atm. 
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VCH}OH  =  0.04072l\liter/mol];  VH0  =0.018068 [liter/mol]. 


From  Eq.  (9 If) 


(l)(0. 018068) 


cH.m^ppr  1  +  (^(0.018068  -  0.040727) 

Difference: 


=  0.01849 


^  Cll -OI  f , exact  ^  CH-X)! !  .appr 


v  CH  3 OH ,  exact 


x  100  =  0.33%  (i.e.  less  than  0.5  %) 


So,  for  practical  purpose,  one  may  use  Eq.  (9 If)  to  compute  mole  fraction  of  CH3OH  in  solution 
instead  of  the  exact  procedure,  requiring  more  work. 


4.2  Plotted  Data 

The  equilibrium  cell  voltage  as  a  function  of  cell  temperature,  pressure,  and  reactant 
concentration  is  given  Figures  2-85  and  starts  on  the  next  page. 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  2.  Set  1  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  3.  Set  1  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  4.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  5.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  6.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  7.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  8.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  9.  Set  1  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  10.  Set  1  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 


Figure  11.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  12.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  13.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  14.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  15.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  16.  Set  1  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  17.  Set  1  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  18.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  19.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


1.200 
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a) 
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ro 
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Figure  20.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  21.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  22.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


temperature  =  25  °C 
anode  pressure  =  1  bar 


c(MeOH)  =  0.01  M 
■«— c(MeOH)  =  0.10  M 
■*—  c(MeOH)  =  0.50  M 


Figure  23.  Set  1  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  24.  Set  1  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  25.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  26.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  27.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  28.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  29.  Set  1  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  1  -  Concentration  of  methanol  range  from  0.01  to  0.50  M. 


Figure  30.  Set  1  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 


Figure  31.  Set  1  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  32.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Figure  33.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  34.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Figure  35.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  36.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  37.  Set  2  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  38.  Set  2  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 


Figure  39.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  40.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  41.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  42.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  43.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


62 


Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  44.  Set  2  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  45.  Set  2  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  46.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  47.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  48.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  49.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  50.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Figure  51.  Set  2  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  52.  Set  2  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  53.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  54.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  55.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  2  -  Concentration  of  methanol  range  from  0.50  to  1.00  M 


Figure  56.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  57.  Set  2  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  58.  Set  3  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 


Figure  59.  Set  3  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  60.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  61.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  62.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Figure  63.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  64.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  65.  Set  3  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  66.  Set  3  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  liquid  phase). 


Figure  67.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  68.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  69.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 
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Figure  70.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 


Figure  71.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  liquid  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  72.  Set  3  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  73.  Set  3  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  74.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  75.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  76.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  77.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  78.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  ideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Figure  79.  Set  3  -  Effect  of  cathode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  80.  Set  3  -  Effect  of  anode  pressure  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas 
behavior  (water  formation  at  cathode  in  vapor  phase). 


Figure  81.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  82.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  83.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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Set  3  -  Concentration  of  methanol  range  from  1.00  to  2.00  M 


Figure  84.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 


Figure  85.  Set  3  -  Effect  of  temperature  on  the  open  circuit  voltage  for  the  nonideal  solution  and  gas  behavior 
(water  formation  at  cathode  in  vapor  phase). 
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5  Conclusions 


The  open-circuit  voltage  equation  for  a  direct  methanol  fuel  cell  was  developed  to  determine  the 
maximum  reversible  cell  voltage  obtainable  as  a  function  of  cell  temperature,  pressure,  and 
methanol  feed  concentration.  Based  on  the  developed  equations  and  generated  figures,  the 
optimal  conditions,  based  on  thermodynamics,  occur  when  the  cathode  pressure  is  high,  the 
anode  pressure  is  low,  the  concentration  of  methanol  is  high,  and  the  cell  temperature  is  low. 
Other  factors  affect  the  optimal  actual  fuel  cell  conditions,  but,  the  effort  here  delineates  the  true 
maximal  thermodynamic  potential  available. 
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